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Abstract 

The high accuracy of modern space astrometry requires the use of General Relativity to model 
the propagation of stellar light through the gravitational field encountered from a source to a given 
observer inside the Solar System. In this sense relativistic astrometry is part of fundamental physics. 
The general relativistic definition of astrometric measurement needs an appropriate use of the 
concept of reference frame, which should then be linked to the conventions of the IAU Resolutions 
[H, which fix the celestial coordinate system. A consistent definition of the astrometric observables 
in the context of General Relativity is also essential to find uniquely the stellar coordinates and 
proper motion, this being the main physical task of the inverse ray tracing problem. Aim of this 
work is to set the level of reciprocal consistency of two relativistic models, GREM and RAMOD 
(Gaia, ESA mission), in order to garantee a physically correct definition of light direction to a star, 
an essential item for deducing the star coordinates and proper motion within the same level of 
measurement accuracy. 
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I. INTRODUCTION 



The correct definition of a physical measurement requires the identification of an appro- 
priate frame of reference. This applies also to the case of the determination of position 
and motion of a star from astrometric observations made from within our Solar System. 



Moreover, modern instruments homed into space-borne astrometric probes like Gaia 



2] and 



SIM [2]] are targeting accuracy at the micro-arsecond level, or higher, thus requiring any as- 
trometric measurement be modelled in a way that light propagation and detection are both 
conceived in a general relativistic framework. One needs, in fact, to solve the relativistic 
equations of the null geodesic which describes the trajectory of a photon emitted by a star 
and detected by an observer with an assigned state of motion. The whole process takes 
place in a geometrical environment generated by an N-body distribution as could be that 
of our Solar System. Essential to the solution of the above astrometric problem, namely an 
inverse ray tracing from observational data, is the identification, as boundary conditions, of 

nnri 

the local observer's line-of-sight defined in a suitable reference frame (see, e.g. |4|, |5|, l6j). 

Summarizing from the references quoted above, the astrometric problem consists in the 
determination, from a prescribed set of observational data (hereafter observables) of the 
astrometric parameters of a star namely its coordinates, parallax, and proper motion. How- 
ever, while in classical (non relativistic) astrometry these quantities are well defined, in 
General Relativity (GR) they must be interpreted consistently with the relativistic frame- 
work of the model. Similarly, the parameters describing the attitude and the center-of-mass 
motion of the satellite need to be defined consistently with the chosen relativistic model. 

At present, three conceptual frameworks are able to treat the astrometric problem at the 
micro-arcsecond level within a relativistic context. 

The first model, named GREM (Gaia Relativsitic Model) and described in Klioner [3], is 
an extension of a seminal study Klioner and Kopeikin [8] conducted in the framework of the 
post-Newtonian (pN) approximation of GR. This model has been formulated according to a 
Parametrized Post Newtonian (PPN) scheme accurate to 1 micro-arcsecond. In this model 
finite dimensions and angular momentum of the bodies of the Solar System are included and 
linked to the motion of the observer in order to consider the effects of parallax, aberration, 
and proper motion. This model is considered as baseline for the Gaia data reduction Jj|. 
The boundary conditions are fixed by the coordinate position of the satellite and imposing 
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the value of c to the modulus of the light direction at past null infinity. The light path 
is solved using a matching technique which links the perturbed internal solution inside the 
near-zone of the Solar System with the (assumed) flat external one. 

Conceptually similar to the above model is the one developed in Kopeikin and Schafer 
[lo| . Using the post-Minkowskian (pM) approximation, Einstein's equations are solved in 
the linear regime expressing the perturbated part of the metric tensor in terms of retarded 
Lienard-Weichert potentials. Later, Kopeikin and Mashhoon included all the relativis- 
ts effects related to the gravitomagnetic field produced by the traslational velocity/spin- 
depedent metric terms. 

Both works, in the pN and pM aproaches, rewrite the null geodesic as function of two 
independent parameters and solve the light trajectory as a straight line (Euclidean geome- 
try) plus integrals, containing the perturbations encountered, from a gravitating source at 
an arbitrary distance from an observer located within the Solar System. This allows one 
to transform the observed light ray in a suitable coordinate direction and to read-off the 
aberrational terms and light deflections effects, evaluated at the point of observation. The 
main difference between the two approximations appears in the computation of the light 
deflection contributions: in the pN scheme by the technique of asymptotic matching, while 
in the pM one by a semi-analytical integration of the equation of light propagation from the 
observer to the source with retarted time as argument. 

The third and last model, RAMOD, is an astrometric model conceived to solve the inverse 
ray-tracing problem in a general relativistic framework not constrained by a priori approx- 



imations 



121 ] . It exploits the concept of a curved geometry as a common background to 



all steps of its functioning and can be extended to whatever accuracy and physical require- 
ments |l3]. RAMOD therefore is not a just a pN model, contrary to how was referenced in 
Moreover, the same parametrization of the pN/pM approximations can be obtained 
in RAMOD if we limit the model accuracy to the milli-arcsecond level 15|. The full devel- 
opment to the micro-arcsecond level imposes to include the h 0i metric terms and to take 
properly into account the retarded distance effects due to the motion of the bodies of the So- 

n 

lar System [5|. At present, the RAMOD full solution requires the numerical integration of a 
set of coupled non linear differential equations (also called "master equations") which allows 
to trace back the light trajectory to the star initial position and which naturally includes all 
the effects due to the curvature of the background geometry. A solution of this system of 
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differential equations contains all the relativistic perturbations suffered by the photon along 
its trajectories due to the intervening gravitational fields. The boundary conditions fixed 
by the astrometric observable as function of an analytical fully relativistic description of the 
satellite allows a unique solution for a stellar position and motion [4J|. 

The first two models, namely the pN and pM ones, though different, take advantage of a 
similar "language" that facilitate their comparison. RAMOD, on the contrary, is formulated 
in a completely different way. This makes its comparison with the former two a difficult 
task. However, since they are used for the Gaia data reduction with the purpose to create 
a catalog of absolute positions and proper motions, any inconsistency in the relativistic 
model(s) would invalidate the quality and reliability of the estimates. This alone is sufficient 
reason for making a theoretical comparison of the two approaches a necessity. 

In this paper we present the first theoretical comparison, showing how it is possible to 
"extract" the aberration terms from the RAMOD construct. 

In section [Til we review all the building steps of the RAMOD astrometric set-up. In 
section Hill we compare the procedures used in GREM to those utilized in RAMOD to define 
the observables and suggest a possible way to make a comparison between the quantities of 
these two formulations via the explicitation of the aberration part in the RAMOD framework. 
Section [TV] is devoted to describe the GREM calculations of stellar aberration, while the 
following one shows how the same effect can be recovered in RAMOD. Section |VT] will 
finally comment on the results of the comparison and on some crucial points which have to 
be addressed to proceed further with the theoretical comparison of the two models. 



II. THE RAMOD FRAMES 



The set-up of any astrometric model implies, primarly, the identification of the gravi- 
tational sources and of the background geometry. Then one needs to label the space-time 
points with a coordinate system. The above steps allow us to fix a reference frame with 
respect to which one describes the light trajectory, the motion of the stars and that of the 
observer. 

The RAMOD framework is based on the weak-field requirement for the background ge- 
ometry, which in turn has to be specialized to the particular case one wants to model. For 
example, having in mind a Gaia-like mission, we can assume the Solar System as the only 
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source of gravity, i.e. a physical system gravitationally bound and weakly relativistic. Then, 
only first order terms in the metric perturbation h (or equivalently in the constant G as in 
the post-Minkowskian approximation) are retained. These terms already include all of the 
possible (t>/c) "-order expansions of post-Newtonian approach, but just those up to (v/c) 3 are 
needed to reach the micro-arcsecond accuracy required for the next generation astrometric 
missions, like e.g. Gaia and SIM. 

With these assumptions the background geometry is given by the following line element 

ds 2 = g af3 dx a dx^ = (r] a p + h aP + O (h 2 )) dx a dx?, 

where 0{h 2 ) collects all non linear terms in h, the coordinates are x° = t, x l = x, x 2 = 
y, x 3 = z, the origin being fixed at the barycenter of the Solar System, and r] a/3 is the 
Minkowskian metric. 

For this reason, any comparison between RAMOD and GREM requires that both use the 
same metric. In the small curvature limit the metric components used in RAMOD are [li| 

£00 = -1 + /*oo + 0(4) g i = h 0l + O(5) gij = 1 + hooSij + O (4) , (1) 

(2) (3) (2) 

where h 00 = 2U/c 2 , h Qi = U l /c 3 , and U and U l are, respectively, the gravitational potential 

(2) (3) 

and the vector potential generated by all the sources inside the Solar System that can be 
chosen according the IAU resolution B1.3 [lj]. The metric of Eq. ([]]) is also adopted by 
GREM. Finally the subscripts indicate the order of (v/c) (e.g. h 0i ~ 0(3) and 0{n) = 

(3) 

0[(v/cY\). 

A. The BCRS 

In the near zone of the Solar System and with the metric (JTJ) , IAU resolutions provide 
the definition of the Barycentric Celestial Reference System (BCRS), and of the Satellite 
Reference System (SRS) These resolutions, as remarked above, are based on the pN 
approximation of GR which is still consistent with RAMOD, since the perturbation h a /3 to 
the Minkowskian metric in (pQ) can be calculated at any desired order of approximations in 
(v/c) inside the Solar System. 

n 

In RAMOD (see de Felice et al. [12|) a BCRS is identified requiring that a smooth family 
of space-like hypersurfaces exists with equation t(x,y,z) = const. The function t can be 



taken as a time coordinate. On each of these t (x, y, z) = const hypersurfaces one can choose 
a set of Cartesian-like coordinates centered at the barycenter of the Solar System (B) and 
running smoothly as parameters along curves which point to distant cosmic sources. The 
latters are chosen to assure that the system is kinematically non-rotating. The parameters 
x, y, z, together with the time coordinate t, provides a basic coordinate representation of 
the space-time. 

Any tensorial quantity will be expressed in terms of coordinate components relative to 
coordinate bases induced by the BCRS. 



B. The local BCRS 



In RAMOD, at any space-time point there exists a unitary four-vector u a which is tangent 
to the world line of a physical observer at rest with respect to the spatial grid of the BCRS 
defined as: 

u a = (-9o y 1/2 5% ={l + ^d t + . (2) 

The totality of these four- vectors over the space-time forms a vector field which is propor- 



tional to a time-like and asymptotically Killing vector field [12j. The proper time measured 
by each of these observers is proportional to the BCRS coordinate time according to equa- 
tion ((21). To the order of accuracy required for Gaia, the rest space of u can be locally 
identified by a spatial triad of unitary and orthogonal vectors whose choice however can 
only be dictated by specific requirements. A natural choice is that of pointing to the local 
coordinate directions chosen of the BCRS (figured]). 

This frame will be called local BCRS] obviously, the local proper time varies as a function 
of the gravitational potential at the observer's position, as can be deduced from equation 
((2J). In the RAMOD formalism this local BCRS is represented by a tetrad whose spatial 
axes (the triad) coincide with the local coordinate axes, but whose origin is the barycenter 
of the satellite. At the O (h 2 ), this triad is |4j 

K = hoaSS + (l " ^) E (3) 

for a = 1, 2, 3. 

In RAMOD any physical measurement refers to the local BCRS. 
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t (x, y, z) = const 




Figure 1: The local observer world line with respect to the BCRS coordinate system. The spatial 
axes of the BCRS point toward distant sources. The dashed lines are the world lines of the observers 
at rest with respect to the barycenter (B). The rest-space (green area) of u locally deviates from the 
space-like hypersurface with equation t[x, y, z) = const by terms of the order of a micro-arcsecond. 

C. The proper reference frame for the satellite 

The proper reference frame of a satellite consists of its rest-space and a clock which 
measures the satellite proper time. 

The tensorial quantity which expresses a proper reference frame of a given observer is 
a tetrad adapted to that observer, namely a set of four unitary mutually orthogonal four- 
vectors Aq, one of which, i.e. Ag, is the observer's four- velocity while the other A„s form a 
spatial triad of space-like four-vectors. Mathematically the tetrad is found as a solution of 
the following system fl^ |: 

Vab = g^KK ( 4 ) 

which allows one to interpret a tetrad frame also as an instantaneous inertial reference 
frame. The solution of pj, always computed w.r.t. the BCRS, is not trivial since it depends 
on the metric at each space-time point along the world line of the observer. The physical 
measurements made by the observer (satellite) represented by such a tetrad are obtained by 
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projecting the appropriate tensorial quantities on the tetrad axes. 

The same measurements can also be defined by splitting the space-time into two subspaces. 
A time-like observer u' a carrying its laboratory is usually represented as a world tube; in 
the case of a non-extended body, the world tube can be restricted to a world line tracing the 
history of the observer's barycenter in the given space-time. At any point P along the world 
line of u' a , and within a sufficiently small neighborhood, it is possible to split the space-time 
into a one-dimensional space and a three-dimensional one 13fl, each space being endowed 
with its own metric, respectively U a p {u') = —u' a u'p and P a p (V) = g a p + u' a u'p. Clearly, 

9aP = U a p (u') + P a p (lb') . (5) 

The space with metric P a p (u r ) is generated by lines which stem orthogonally to the world- 
line of u' at P and is denoted as the rest-space of the observer u' at P. In this space one 
measures proper lengths. The space with metric U a p (u f ) is generated by lines which differ 
from that of u' by a riparametrization. In this space one measures the observer's proper 
time. 

As a consequence of Equation (jSJ), the invariant interval between two events in space- 
time can be written as ds 2 = P a/ 3 (u 1 ) dx a dx /3 + U a p (V) dx a dx 13 , from which we are able to 
extract the measurements of infinitesimal spatial distances and times intervals taken by u' 
as, respectively, dL u > = \JP a p (V) dx a dx 13 and 

dT u , = -c- l u' a dx a . (6) 

Essentially, the last method is equivalent to the tetrad formalism, when we do not know 
the solution of (j4j) and we need to know only the moduli of the physical quantities. As 
far as RAMOD is concerned, given the metric (pQ) and in the case of a Gaia-like mission, 
an explicit analitic expression for a tetrad adapted to the satellite four-velocity exists and 
can be found in [J]. The spatial axes of this tetrad are used to model the attitude of the 
satellite. Moreover, from eq. ((6l), it is possibile to deduce the IAU trasformations between 
the observer's proper time and the barycentric coordinate time, without using any matching 
tecnique This finally sets the running time on board and completes the definition of 

the proper reference frame for the Gaia-like satellite. 
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III. SINGLE-STEP VS MULTI-STEP DEFINITION OF THE OBSERVABLE AND 
A WAY FOR THE RAMOD VS. GREM COMPARISON 



The classical (non relativistic) approach of astrometry has traditionally privileged a 
"multi-step" definition of the observable; i.e., the quantities which ultimately enter the "fi- 
nal" catalogue and are referred to a global inertial reference system, are obtained taking 
into account, one by one and independently from each other, effects such as aberration and 
parallax. 

GREM reproduces in a relativistic framework this approach of classical astrometry. The 
BCRS is, for this model, the equivalent of the inertial reference system of the classical 
approach, while the final expression of the star direction in the BCRS is obtained after 
converting the observed direction into coordinate ones in several steps which divide the 

n 

effects of the aberration, the gravitational deflection, the parallax, and proper motion f7j. 

In the previous section we have mentioned that RAMOD relies on the tetrad formalism 
for the definition of the observable. In general, the three direction cosines which identify 
the local line-of-sight to the observed object are relative to a spatial triad associated to 
a given observer u'; the direction cosines w.r.t. the axes of this triad are defined as: 

P{u') a pk a El 

cosipa = %-^ = e & , (7) 

where k a is the four-vector tangent to the null geodesic connecting the star to the observer, 
and all the quantities are obviously computed at the event of the observation. 

As a consequence of this definition, given the solution of the null geodesic equation and 
the motion and the attitude of the observer, equation ((7]) expresses a relation between the 
unknowns, namely the position and motion of the star, and the observable quantities which 
includes all of the above effects mentioned for GREM. In other words, in RAMOD it is 
not needed and not natural to disentangle each single effect, relativistic or not. For this 
reason any attempt to make a theoretical comparison between the two models is difficult, 
but the way how the observer tetrad was found in RAMOD suggests a way to overcome this 
problem. 

In Bini et al. [J| the attitude frame E~ a was strictly specified for measurements made by 
a Gaia-like observer. Let us summarize the main steps. 

Given the tetrad {A^} adapted to the local barycentric observer as defined in (Bini et al. 
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4] and reference therein) the vectors of the triad {Aa} are boosted to the satellite rest frame 
by means of an instantaneous Lorentz transformation which depends on the relative spatial 
velocity v a of the satellite identified by the four-velocity u' a w.r.t. the local BCRS u a , and 



whose Lorentz factor is given by 7 = —u a u' a [17]. 



The boosted tetrad A a obtained in this way represents, simarly to what is defined for 
Gaia in Q], a CoMRS (Center-of-Mass Reference System, comoving with the satellite). 
In addition to the definition in the cited works one of the axes is Sun-locked, i. e. one axis 
points toward the Sun at any point of its Lissajous orbit around L2, in order to deduce the 
Gaia attitude frame. This final task is obtained by applying the following rotations to the 
Sun-locked frame: 

1. by an angle uj p t about the vector A" which points constantly towards the Sun, where 

bs J 

uj p is the angular velocity of precession; 

2. by a fixed angle a about the image of the vector A" after the previous rotation; 

bs 2 " 

3. by an angle u r t about the image of the vector A" after the previous two rotations, 

bs ] 

where u r is now the spin angular velocity. 
The triad resulting from these three steps establishes the satellite attitude triad, given by: 

E h = TZ l (u r t)TZ 2 (a)TZ 1 (uj p t)\ a = 1,2, 3. 

bsa 

The final triad should be the equivalent, in the RAMOD formalism, to the Satellite 
Reference System (SRS) £jj of GREM. 

Once this procedure is completed, the final measurements will naturally entangle in a 
single result every G R "effect". Therefore, the natural way to "extract" any of those effects 
in a separate formula, is to consider equation (jjj) and express the observable as a function 
of the appropriate tetrad. 

IV. STELLAR ABERRATION IN GREM 

As well known stellar aberration arises from the motion of the observer relative to the 
BCRS origin, assumed to coincide with the center of mass of the Solar System. In order to 
account for stellar aberration in the algorithm for the reduction of the astrometric observa- 



tions, the pN/pM approaches 
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ll| transform the observed direction to the source (s) into 
10 



the BCRS spatial coordinate direction of the light ray at the point of observation x(£ Q ) = x s 
(see figure [2j). Now, paraphrasing Klioner 2J, the coordinate direction to the light source at 
x s is defined by the four- vector p a = where p % = c~ x dx % /dt, x % and t being the BCRS 

coordinates. But the coordinate components p l are not a directly observable quantities; the 
observed vector towards the light source is the four- vector s a = (1, defined with respect 
to the local inertial frame of the observer. In the local frame: 

(8) 



where X a are the coordinates in the CoMRS, then in order to deduce the spatial direction 
p l from s l it is chosen to proceed as follows. 

From the property of a null trajectory and taking into account the metric which defines 
the BCRS it is 



namely 



26 {j (-2(1 +7ppn) 



-1 + 



(1 + 7p PN ) w(t, x) 2w 2 (t,x) 



w l (t, x) 



V 3 + $ij 1 + 



;i + 7ppn)^(^,^) 



+ 

pV = 



which gives 



1 _ (1 +7ppn) w(t,x) _ 2 (l+'jp PN )5 ij w l (t,x)p j _ 4 
p C A c 6 



(9) 



where p = a/ 5ijp l pi is the Euclidean modulus of the spatial vector p and 7ppn is the PPN 
parameter. 

The infinitesimal transformation X a (x^) between CoMRS and BCRS is given by the 
formula: 

dX a = A^daA (10) 
From (fTUl ) the expression of s l as a function of the spatial components p % is obtained: 



Ag + AjpJ" 



(11) 
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star path 




Figure 2: The vectors representing the light direction in the pM/pN approches inside the near-zone 
of the solar system. 

One can explicit formula (TTT1) by following the procedure reported in [8| and adopting the 
IAU resolution B1.3 l\. From the BCRS (ct, x l ) to the CoMRS (cT, Af*) s the transformation 
between the time coordinates reads: 



T = t - c~ 2 [A(t) + s^Ri] 

+c- 4 [B + SijB'Rl + 5 im 5 ]k B %1 R™R k s + C(t, x)] + O (c~ 5 ) , 



(12) 
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and between the spatial coordinates 



5) + c~ 2 Q^ + ff^W +*>}(*)) 



Hi ■ <■ Vr^niinit ■ o(<- ; ) . (13) 



All the functions A, B, C, D are defined in Klioner and Kopeikin [8| or in IAU resolutions 
and 

are the coordinate displacements with respect to the center of mass of the satellite x l s in the 
BCRS, and finally 

dt 

is the coordinate velocity of the center of mass of the satellite relative to the BCRS. 

As reported in Klioner [lj3], the attitude in GREM (SRS) is obtained by applying an 
orthogonal rotation matrix 71° to X 1 in equation f fT3"l) . At this stage the role of the SRS is 
equivalent to that of the E?s in eq. (jZj). 

If one keeps all the terms up to the order of 1 micro-arcsecond, the observed coordinate 
direction s l , in terms of the unitary spatial vector n l = p l /p, becomes in the CoMRS: 

s l = —n % + cT 1 [n x (v x n)f 

+cT 2 |(n ■ v) [n x (n x v)f + ^ [v x (n x v)f j 

+c" 3 | [(v • n) 2 + (1 + 7 ppnV (x s )] [n x (v x n)] 1 



+i(n.v)[vx(nxv)] i }+0&. (14) 



V. RAMOD ABERRATION IN THE PM APPROXIMATION 

Whatever tetrad we consider, the expression of Eq. j7|) for the relativistic observable in 
the RAMOD model can also be written as Q| 



7 (l " vJfo)) 



(15) 



where v a is the spatial four- velocity (also called as the "physical velocity") of the satellite u' 
relative to the local baricentric observer u. The quantity l^ was introduced in RAMOD 
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and is a unitary four-vector which represents the local line-of-sight of the photon as seen by 



u, i.e. I 



(o) 



PSitfkP. 



Finally, 7 is the Lorentz factor of u' a with respect to u a , that is, 



— u' a u a 



1 



= 7, 



(16) 



where v 2 = v a v a . 

To retrieve the aberration effect given by the motion of the satellite with respect to the 



BCRS in RAMOD, one needs to specialize Eq. ( jT5h to the case of a tetrad |Aq j adapted to 
the center of mass of the satellite assumed with no attitude parameters. In this case, in fact, 
the observation equation will give a relation between the "aberrated" direction represented by 
the direction cosines cosipa as measured by the satellite and the "aberration-free" direction 
given by the quantity l? s referred to the local BCRS frame A?. The vectors of the triad 
Xa \ differ from the local BCRS's {Aa} for a boost transformation with four-velocity u' a . 



This means that it can be derived from Eq. ((3j using the relation 
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X? = P(u' 



A? 



7 



7 + 1 



v a KA p a) 



(17) 



where u' a and u a are the above mentioned four-velocity of the satellite and its physical 
velocity relative to the local BCRS respectively, and P («')" = + u' a u' a . 
From de Felice et al. Q| and Bini et al. [l] it is 



u 



(18) 



where v l is the coordinate velocity of the satellite, as stated in the previous section. Now, 
being |j| 

v a = - (u' a - 1U a ) , (19) 
7 

one deduces that u° ~ O (v 4 / c 4 ) and 



U\ v l _ / v 4 
+ O 



v 

cf J c yc" 1 
Expanding Eq. ( ITTjl with relations ( TT61) and (fT9l one gets 



(20) 



Aa = A? + u' a («'|A a ) 



(21) 
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where the notation ( 
Then, using Eqs. 
we obtain 



represents the scalar product, so, e.g., (w'|Aa) = g a pu' a \^. 

lfl8~l) and (l20l) and expanding the scalar products to the right order 



U'\\a) 


-0- 


(^|A a ) 


-0 






(u'\u) 





2t/ 1 v 2 \ v a 



2 c 2 / c 



+ o 



2U\ v 



so that the expression for the boosted tetrad finally becomes 

A" - A" + (l + 3f/ + 1 8^ + 1 "V*" + O ^ 4 



(22) 
(23) 
(24) 

(25) 



A. The (v/c) expansion of the relativistic observable 



Given Eq. (1251) one can consistently recast Eq. (11511 as 



p 

13 °0 



^ 7 (1 - ' 7 (1 - V c 2 + 2 c 2 ) c + 2 7 (1 - vj°) c W ' 



where e a are the cosines related to the tetrad A? which, as said, does not contain the attitude 
parameters. Here and in the rest of the section, we replace the symbol l? \ with l a to ease 
the notation. 

After long calculations, the first term on the right-hand-side of this formula can be written 

as 



r + - \-v a + (SijvH*) r] + 

7 (1 - uj a ) c 1 v ' J 



ur - [Sijv'v) v a + 



7i\2 



— <^ -2f/w a - 



2 1 



r + 



3Z7 (>V'7 ? ) + [SijvH^Y - ^v 2 (S^vV) + U (S^P) 



the second term is zero since both /oand ^ are zero, while the third one becomes 



1 (l-iy).(v i /c) v a 1 

2 7 (1 - ~ = 2 



Sid 1 h Stf—l 3 



. 2 o 
!)'-.\ f 



V IV 

— + o - 
c V c 



(27) 



(28) 



15 



Finally, collecting all terms: 



e a = l a + ~ [~ va + {Sijv'F) l a ] + 
lSui a -l(S ij v i F)v a + 



2 1 2 
V 

2 



r } + 



^{-2C/^-I(^V)V+ 



r 



3U (SijvHi) + (Sij-vH*) - -v 2 (SijvH*) + U (S^P) 



O 



(29) 



At a first glance, the last expression shows differences in terms up to the (v/c) 2 order (note 
in particular the appearance of the term Ul a ) and of the (v/c) 3 order which cannot allow to 
straightforwardly compare, as expected, the above expression to the GREM vectorial one of 
Eq. (HD. 



B. Comparison with the GREM model 

The expression (1291) relates the observed direction cosines with l a . The equivalent relation 
for the GREM observable is equation (fT4l) where the aberration is expressed in terms of a 
vector n. To compare formula (|29l) with GREM's formula (IT3|) we need to find a relationship 
between n and l a . To this purpose we need to reduce T a to its coordinate euclidean expression 

In GREM n represents the "aberration-free" coordinate line of sight of the observed star 
at the position of the satellite momentarily at rest. In RAMOD, as said, l a represents the 
normalized local line- of- sight of the observed star as seen by the local barycentric observer 
u. In other words, l a is a four-vector which fixes the line-of-sight of an object with respect 
to the local BCRS. 

Do n and l a have a similar role in the two approaches? From the physical point of view 
they have the same meaning, as the observed "aberration free" direction to the star. Let us 
start from the definition of n in GREM: 

i P i 

71=—, 

V 

where p l = c _1 dx l /dt and p is the Euclidean norm of p\ so that p^ 1 ~ (1 + h 00 + h 0i p l ) + 
O (/i 2 ), as equation ([9l) shows. This means that 

n l =p i (l + /ioo + W) + O (h 2 ) . (30) 
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On the other hand, using the definition of l a in de Felice et al. [eJ it can be easily shown 
that its spatial components are 

ji = = i 

1 ^ 2 and k l /k° = c _1 dx*/dt = p\ it results 
1/2 /, , , 



and, from u = (—goo 

T* = P l (-9o ) 1/2 (1 - h 00 - W) 
= p* (l + hoo + Wj + O (h 2 ) . 

Finally, from equations ( 1301) and (I3TT) one has 

f = n* ( 1 - h I O 



(31) 



(32) 



namely, the spatial light direction, expressed in terms of its Euclidean counterpart at the 
satellite location in the gravitational field of the solar system. Worth noticing is that no 
terms of the order of 0[(v/c) 3 } appear in (l32l) . 

Combining Eq. (l29l) with (l32l) and setting (SyVW) = v • n to ease the notation, we 
obtained 



1 1 
n a + - [-v a + (v • n) n a ] + — 



1 



v ■ n) v a + 



v • n 



n 



+ 



4 < -2Ur" 



-(vn) 2 !)' 
2 v ' 



v • n n 



2C/ + (v-n) 



1 s 

-IT 

2 



H-O -T ■ (33) 



In this way the right-hand side of the aberration expression of RAM OD is rewritten with 
the GREM quantities at the (f/c) 3 order. The same operation can be done for the left-hand 
side using the definition of the projection operator and the tetrad property A^A^g = g a p: 

P(u) a(3 k a \ p , k a Xt k a \i k a \1 dx a 



(34) 



Is there a relation between the direction cosines of the above equation with the spatial 
components of the observed vector s l in GREM? The crucial point stands on the definition of 
the coordinates system. The tetrad components of the light ray can be directly associated to 
CoMRS coordinates (as done in Klioner [l^|) if the boosted local BCRS tetrad coordinates 
x a are equivalent to the CoMRS ones X a . This is true only locally, i.e. in a sufficiently 
small neighborhood (since the tetrad are not in general olonomous) and if the origins of the 
two reference systems concide. So, from (JSj) , if one could state that 
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dx & _ dX a 
~ dX° 



it would follow 



(35) 



In RAMOD, at the milli-arcsecond level, the rest space of the local baricentric observer 
coincides globally with the spatial hypersurfaces which foliate the space-time and define 



the BCRS 



Y2\. At micro-arcsecond accuracy, instead, the vorticity cannot be neglected 



and the geometry is affected by non-diagonal terms of the metric hence the t = constant 
hypersurfaces do not coincide with the rest-space of the local barycentric observer p|. Then, 
to be consistent we can only define at each point of observation a spatial direction measured 
by the local barycentric observer and then associate it to the satellite measurements via the 
direction cosines relative to the boosted attitude frame. As far as GREM is concerned, the 
euclidean geometry admits a parallel transport which does not feel the curvature, allowing 
to define the same vector in any point of the space. 

Then, equation f[35T > has only local validity and (|33l)can be written as 



n a + - [-v a + (v ■ n) n a ] + 4 i -- (v • n) r" 



v ■ n 



4 I -2Uv a - - (v • n) 2 v a + (v • n) n a 



2U+ (v-n) s 



—v 
2 



-(36) 



Considering that n • n = 1 and v 2 = 5ijV % iP = v ■ v, the previous equation becomes 



s a = ~n a + - [v a (n ■ n) - n a (v ■ n)] + 



1 ^ (v • n) [v a (n ■ n) - n a (v ■ n)] + ^ [n a (v ■ v) - v a (v ■ n)] | + 

{2U [v a (n • n) - n a (v • n)] + (v • n) 2 \v a (n • n) - n a (v • n)] + 

1 1 (v A 

(v ■ n) [n a (v • v) - v a (v • n)] } + O 



2 , - v .,. (3T) 

Finally, from the relation ax (b x c) = b (a • c) — c (a ■ b) it is 

s a = —n a + - [n x (v x n)] a + < (v • n) [n x (v x n)] a + - [v x (n x v)] a 1 + 

i | [(v • n) 2 + 2U] [n x (v x n)f + \ (v ■ n) [v x (n x v)f} + O (£\ (38) 

which is formula (TT4ll for the aberration in GREM if we consider the case of GR where 
7ppn = 1 and we take into account that v = x D , and U = w (x Q ). 
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Finally, the result obtained with eq. (l38l) states that, limited to the case of aberration and 
using the appropriate definitions of the IAU recommendations, RAMOD recovers GREM at 
the (f/c) 3 order. 



VI. CONCLUSIONS 

This paper compares two relativistic astrometric models, GREM and RAMOD, both 
suitable for modelling modern astrometric observations at the micro-arcsecond accuracy. 
Their different mathematical structures hinder a straightforward comparison and call for a 
more in-depth analysis of the two models. Because of the structure of GREM, the earliest 
stage of a theoretical comparison starts with the evaluation of the aberration "effect" in 
RAMOD. In this regard, we can evidence the following differences in: (i) the choice of the 
boundary conditions, (ii) the tools needed to define the astrometric measurements, (iii) the 
attitude implementation, (iv) the definition of the proper light direction. 

Crucial is point (i). The light signal arriving at the local BCRS along the spatial direction 
l a = Piu^k 13 satisfies the RAMOD master equations, namely a set of non-linear coupled 
differential equations {5]. Therefore the cosines (i.e. the astrometric measurements) taken 
as a function of the local line-of-sight (the physical one), at the time of observation (/( )), 
allow to fix the boundary conditions needed to solve the master equations and to determine 
uniquely the star coordinates. However, since the direction cosines are expressed in terms 
of the attitude, the mathematical characterization of the attitude frame is essential to 
complete the boundary value problem in the process of reconstructing the light trajectory. 
The vector n, i.e. the "aberration-free" counterpart of l a in GREM, is instead used to derive 
the aberration effect (in a coordinate language) and there is no need to connect it with a 
RAMOD-like boundary value problem. 

As for the solution of the geodesic equation, RAMOD defines a complete procedure to 
derive the satellite attitude which depends as input only on the specific terms of the metric 
that describes the addressed physical problem. GREM, instead, embeds the definitions of its 
main reference system (BCRS) within the metric, consequently each further step depends on 
this choice. This includes all the subsequents transformations among the reference systems 
which are essential to extract the GREM observable as function of the astrometric unknowns. 
On the other side, the RAMOD analytical solution for the attitude frame assures controlled 
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alghoritms that can be directly implemented in the solution of the astrometric problem and 
guarantee its consistency with GR. In RAMOD the direction cosines link the attitude of the 
satellite to the measurements, compacting several reference frames useful to determine, as 
final task, the stellar coordinates: the BCRS (kinematically non-rotating global reference 
rame), the CoMRS (a local reference frame comoving with the satellite centre of mass), 
and the SRS (the attitude triad of the satellite). The coordinate transformations between 
BCRS/CoMRS/SRS come out naturally once the IAU conventions are adopted. This is 
inside the conceptual framework of RAMOD, where the astrometric set-up allows to trace 
back the light ray to the emitting star in a curved geometry, and it is not natural to disen- 
tangle each single effect. Any approximation can be applied a posteriori where it is needed, 
case by case. This explains items (ii) and (iii) and introduces item (iv). 

The direction cosines being physical quantities not depending on the coordinates, are a 
powerful tool to compare the astrometric relativistic models: their physical meaning allow 
us to correctly intepret the astrometric parameters in terms of coordinate quantities. This 
justified the conversion of the physical stellar proper direction of RAMOD into its analgous 
Euclidean coordinate counterpart, which ultimately leads to the derivation of a GREM- 
style aberration formula. Another point arises when the observables of RAMOD have to be 
identified with components of the observed s l of GREM. This matching is admitted only 
if the origins of the boosted local BCRS tetrad in RAMOD and of the CoMRS in GREM 
concide. 

To what extent the process of star coordinate "reconstruction" is consistent with 
GR&Theory of Measurements? Solving the astrometric problem in practice means to com- 
pile an astrometric catalogue at same order of accuracy of the measurements. This paper 
shows that, already at the level of the aberration effect, a correct treatment of physical 
meaurements in terms of coordinate quantities needs particular care in order to avoid mis- 
understandings in the interpretation of the quantities which constitute the final catalogue. 

As a closing consideration, the computation of the BCRS stellar direction in GREM 
needs to extract, at a second stage, the deflection terms from the coordinate "aberration- 
free" direction n l . This problem in RAMOD is, again, embedded in the formulation of the 
astrometric problem as "global solution" which aims at recovering the star coordinates by 
integration of the geodesic equations (treated with an appropriate physical boundary condi- 
tion and approriate reference systems) where the deflection terms play the most fundamental 
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role. 
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